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ABSTRACT 

Recent prope r motion measurements o f the Large and Small Magellanic Clouds (LMC and SMC, 
respectively) by iKallivavalil et al.1 (|2006alTbl ) suggest that the 3D velocities of the Clouds are substan- 
tially higher (~100 km/s) than previously estimated and now approach the escape velocity of the 
Milky Way (MW). Previous studies have also assumed that the Milky Way can be adequately mod- 
eled as an isothermal sphere to large distances. Here we re-examine the orbital history of the Clouds 
using the new veloc ities and a ACDM-motivated MW model with virial mass M V i r = 10 12 M Q (e.g. 
iKlvoin et all (l2003 )L We conclude that the L/SMC are either currently on their first passage about 
the MW or, if the MW can be accurately modeled by an isothermal sphere to distances > 200 kpc 
(i.e., M V i r > 2 x 10 12 M Q ), that their orbital period and apogalacticon distance must be a factor of 
two larger than previously estimated, increasing to 3 Gyr and 200 kpc, respectively. A first passage 
scenario is consistent with the fact that the LMC and SMC appear to be outliers when compared to 
other satellite galaxies of the MW: they are irregular in appearance and are moving faster. We discuss 
the implications of this orbital analysis for our understanding of the star formation history, the nature 
of the warp in the MW disk and the origin of the Magellanic Stream (MS), a band of HI gas trailing 
the LMC and SMC that extends ^100 degrees across the sky. Specifically, as a consequence of the 
new orbital history of the Clouds, the origin of the MS may not be explainable by current tidal and 
ram pressure stripping models. 

Subject headings: galaxies: interactions — galaxies: kinematics and dynamics — galaxies: evolution 
- Galaxy: structure — Magellanic Clouds 



1. INTRODUCTION 

Accurate estimates of the current positions and veloc- 
ities of the Large and Small Magellanic Clouds (LMC 
and SMC, respectively) are crucial parameters for mod- 
eling their past interactions with the Milky Way (MW). 
Recent proper motion measurements with the Hub- 
ble Syace Telescoye (HS T) of the LMC and SMC by 
IKallivavalil et all (|2006allbl . hereafter Kl and K2) im- 
ply a substantial increase (~ 100 km/s) in the estimate 
of the current 3D galactocentric velocities of the Clouds. 
This paper explores the consequences of these new mea- 
surements for our understanding of the global dynamics 
of the MW-LMC-SMC system. 

Constraints on the orbital history of the Clouds can be 
obtained from the structural and kinematic properties of 
the Magellanic Stream (MS). The MS is a band of HI 
gas structured in filaments and clumps extending ^100° 
across the sky, nearly along a great circle beginning at 
the Clouds and passing through the so uth galactic pole 
(|Bruns et al. 2005; Put man et alJl2003T ). From radial ve- 
locity measurements of the HI gas, the Clou ds are known 
to be leading the MS (|Putman et al.lll998t) . Also, since 
the potential well of the Clouds is shallow compared to 



that of the MW, the orbit of the Clouds sho uld be simi- 
lar to that of the MS (jFich fc Tremainelll99lh . Thus, the 
Clouds are believed to be currently following a nearly po- 
lar orbit of counter-clockwise sense as seen from the Sun. 

Previous studies of the orbital evol ution of the Clouds 
and t he formation of the M S (e.g. iMurai fc Fuiimotol 
(|1980l hereafter MF80L |Lin fc Lvnden-Belll (11982T 
LL82L [Gardiner et all (11994 G SF94L Helle r fc Rohlfel 
(fl99l HR94l.lMoore fc Davisl (fl99l. lLin et all (fl995h 
iGardiner fc Noeuchil (I1996LG N96L iBekki fc Chiba 
' 2005h . lYoshizaw aJ^Ngguchj(l2003T ). IConnors et"al 



2005h . iMastropietro et all (|2005l . M05)) generally fol- 
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low the prescription introduced by MF80: the differential 
equations of motion for both the LMC and SMC are in- 
tegrated numerically, allowing the position and velocity 
of the Clouds to be followed backwards in time. Because 
of the large uncertainties in previous proper motion es- 
timates for the Clouds vdM02]vanderMarel, the listed 
authors chose orbital parameters that best reproduced 
the properties of the MS under the assumption that the 
LMC and SMC form a binary system that has been in a 
slowly decaying orbit about the MW for roughly a Hub- 
ble time (£h)- 

The existence of a common envelope of diffuse HI gas 
surrounding the Clouds supports the a ssumption of bina- 
rity fo r at least some time in the past (jMathewson et aU 
1974). However, the only constraint on the orbital pe- 
riod of this binary system about the MW is the length of 
time required to form the MS. This timescale is largely 
dependent on the proposed formation mechanisms: gen- 
erally, theories invoke some combination of ram pressure 
and tidal stripping, implicitly assuming that the Clouds 
have undergone multiple pericentric passages. These ar- 
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guments will be revisited critically in £15.41 

Recently, Kl and K2 have measured the systemic 
proper motions of the LMC and SMC by using the Ad- 
vanced Camera for Surveys (ACS) on HST to track the 
LMC (SMC) 's motion relative to 21(5) background QSOs 
discovered from their optical variability in the MACHO 
database. This has made it possible to determine the 
proper motion of the LMC to better than 5% accuracy: 
[i w = -2.03 ± 0.08mas/yr, fi N = 0.44 ± 0.05mas/yr 
(Kl). They also obtain a four- fold improvement in ac- 
curacy over previous measurements of the SMC's proper 
motion: — —1-16 ± 0.18 mas/yr, //jy = —1.17 ± 0.18 
mas/yr (K2). The corresponding mean 3D space veloci- 
ties and one-sigma errors for the LMC are listed in row 
four of Table [TJ A lower limit on the interaction time 
between the LMC/SMC binary and MW is thus -250 
Myr, which corresponds to the amount of time it takes 
the LMC to travel the length of the MS, assuming that 
the MS extends along a circle of radius 55 kpc and that 
the LMC is moving at v = 378 km/s (the mean 3D ve- 
locity determined by Kl). 

Previous estimates of the LMC's proper motions are 
summarized in Table O The Kl proper motions differ 
substantially from those used in the theoretical mod- 
els of GN96. However, only the west component of the 
Kl values differs markedly from the weighted average of 
previous measurements determined by vdM02. In this 
analysis we focus mainly on the HST measurements, al- 
though we also consider the vdM02 proper motion values 
as a point of reference. We have not explicitly included 
the recent iPedreros et al.l (|2006[ ) results in our analysis 
as this work is essen tially a refinement of the previous 
IPedreros et alJ (|2002l ) measurements, which are already 
included in the vdM02 weighted average. 

The significant error reduction in the HST proper mo- 
tion measurements has also limited the range of plausi- 
ble current orbital parameters. Specifically, the 4a lower 
limit on the LMC's present 3D velocity is now — 310 
km/s, which is larger than the GN96 estimate of 297 
km/s. Correspondingly, if we model the MW as an 
isothermal sphere, the predicted apogalacticon distance 
reached by the LMC on its last orbit would be roughly 
twice as large are previously predicted (see § [2] and Ta- 
ble [TJ). This change owes to the substantial increase in the 
estimate of the LMC's tangential velocity component (i.e. 
Hw)- The large tangential velocity (vtan) and relatively 
small radial velocity (v ra d) imply that the Clouds are 
likely either near apogal acticon or perigalacticon. Ear- 
lier works, such as that o f iLvnden-Bell fc Lin! (|1977f ) and 
iFujimoto fc Sofuel (|1976f ). assumed that the Clouds were 
near apogalacticon; however, the debate was effectively 
settled by early proper motion measurements which un- 
equivocally showed that v tan > 220 km/s (i.e. larger 
than roughly the maximum circular velocity). The cen- 
trifugal force is thus larger than the gravitational force 
acting on the LMC, meaning the LMC is moving toward 
larger radii [dr/dt > 0) and so must be at perigalacticon. 

Yet, a tangential velocity comparable to that measured 
by Kl (367 km/s) has been considered by just three stud- 
ies before: MF80, LL82 and HR94 (v tan = 340, 370 and 
350 km/s, respectively). They determined the orbital pa- 
rameters of the LMC by requiring that its projected orbit 
on the sky traced the current location of the MS (see ^4.2|) 
and that the Clouds have been bound to one another for 



~ 10 10 Gyr. However, GSF94 made the same assump- 
tions and yet determined a significantly lower tangen- 
tial velocity (vtan = 287 km/s). The requirement that 
the LMC's projected orbit traces the MS only uniquely 
determines fijsr] there exists a substantial degeneracy in 
possible nw values when a 3D orbit is viewed in pro- 
jection on the plane of the sky (see £)4.2|) . Since it is 
Hw that ultimately determines the tangential velocity 
component, it becomes relatively clear how these groups 
could make the same assumptions and obtain such differ- 
ent estimates for v tan - For example, although the values 
of /x/v determined by the theoretical models of HR94, 
LL82 and GSF94 are roughly the same (Table EJ), their 
estimates for /j,w are quite different. GSF94 also matched 
the line-of-sight radial velocities along their orbit to the 
available HI data in the MS in determining their proper 
motions; but, this assumption is dependent on the ad- 
vocated formation mechanism of the MS, which is still a 
subject of debate (see §5.4.3j) . The discrepancy between 
the fiw components in these studies ultimately owes to 
the details of the MW models and choice of perigalac- 
ticon, which strongly affects the lifetime of the binary 
system. 

From Table [TJ only the HR94 orbital parameters are 
consistent with those of the Kl orbits in an isothermal 
sphere model. However, as the apogalacticon distance in- 
creases, the validity of an isothermal sphere model for the 
MW is questionable: simulated DM halos don't behave 
like isothermal spheres, especially at large radii. Further- 
more, the corresponding increase in orbital period lim- 
its the effectiveness of ram pressure and tidal stripping, 
which are two popular mechanisms for the formation of 
the MS. 

In this paper we revisit this classic orbital mechan- 
ics problem, adopting 3D velocities constrained by the 
recent proper motion measurements of K1(K2) for the 
LMC(SMC) and a more detailed MW model. In g2]we 
first reconsider the isothermal sphere model to highlight 
the substantial change in the LMC's orbital history im- 
plied by the new measurements. In $3] we consider in- 
stead a cosmologically-motivated model and describe our 
methodology and model parameter s: our fiducial MW 
model is consistent with model Ai of Klv pin et al.l (2002, 
hereafter KZS02) (virial mass M vir = 1O 12 M ) and 
known observational constraints. We focus our analysis 
on the LMC alone as the LMC:SMC mass ratio of -10:1 
precludes the SMC from being a major determinant in 
the LMC's orbital history. In fJTJthe error space of the 
new proper motion measurements is searched to test all 
allowed orbital histories for the LMC. For our fiducial 
MW model, the new mean 3D velocity is roughly the 
escape velocity at the LMC's current location (50 kpc). 
Consequently, our computed orbital histories differ sub- 
stantially from the predictions of all previous studies: we 
find that the orbit of the LMC is best described as highly 
eccentric and approaching parabolic, implying that the 
LMC is on its first passage about the MW. We further 
show that orbits computed with either our fiducial model 
or the isothermal sphere model are co-located when pro- 
jected on the plane of the sky, but deviate from the cur- 
rent location of the MS. Moreover, this deviation exists 
even if the vdM02 proper motion measurements are used. 
In £14.31 our results are shown to be robust to changes in 
our fiducial model parameters. We also consider the im- 
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TABLE 1 

Summary of LMC Orbital Parameters Adopted in Previous Studies 



Work 


3D v (x,y,z) (km/s) a 


\v\ (km/s) 


titan (km/s) 


■u rad (km/s) 


r (x,y,z) (kpc) 


T (Gyr) b 


Peri (kpc) 


Apo (kpc) 


MF80 


(233.7,-13.1, 252. 4) c 


344 


340 


92 


(42.9,-2.4,-28.3) 


1.5 


50 


110 


GSF94, GN96 


(-5,-226, 194) c 


297 


287 


82 


(-1.0,-40.8,-26.8) 


1.5 


45 


120 


HR94 


(-10.06,-287.09,229.73) c 


367.83 


351.81 


107.37 


(-0.85,-40.85,-27.95) 


2.5 


46.3 


180 


vdM02 


(-56 ± 36, -219 ± 23, 


293 ± 39 


281 ± 41 


84 ± 7 


(-0.8, -41.5, -26.9) 


2 


45 


110 




186 ± 35) d 
















M05 e 


(-4.3,-182.45, 169.8) d 


249.3 


237.9 


74.4 


(0,-43.9,-25.04) 


2 


45 


115 


Kl Mean 


(-86 ± 12, -268 ± 11, 


378 ±18 


367 ± 18 


89 ± 4 


(-0.8,-41.5,-26.9) 


3 


50 


220 




252 ± 16) d 
















K2 Fig. 12 


(-91, -250, 220) d 


345 


333 


92 


(-0.8,-41.5,-26.9) 


2 


50 


150 



a All positions and velocities are measured in Galactocentric coordinates. T is the orbital period These velocities are predictions from theoretical 
models.^ These velocities are derived from observations of the LMC's proper motions using the distance moduli of 18.50 ± 0.1 jFreedman et al.H200ll) . la 
errors are quoted. Parameters for M05 were obtained via private communication (2007). 



TABLE 2 

Summary of LMC Proper Motions Estimates 



Work fiyy (mas/yr) /xjv (mas/yr) Method 

GSF94,GN96 -1.72 0.12 MS model 

HR94 a -2.0 0.16 MS model 

vdM02 -1.68 ± 0.16 0.34 ± 0.16 Compilation of pre-2002 measurements 11 

Kl -2.03 ± 0.08 0.44 ± 0.05 HST measurements 



a These are our estimates of the LMC's proper motions corresponding to the 3D veloci- 
ties used by HR94, determined using the distance moduli from IFreedman et al.l d2001f) . b The 
vdM02 value is the weighted average of the ground-based and Hipparcos measurements by: 
Kroupa et al.l U994I) ; IJones et all U994D ; IKroupa fe Bastianl II 19971) ; [Pcdrcros ct al. (2002); 
Drake ct al. (2001). 



pact of the SMC on the robustness of our results in cj5.ll 
We discuss the reliability of the new proper motion mea- 
surements in §5.21 and comment on the likelihood of a 
first passage scenario in £)5.3I The implications of the 
new orbital history for the star formation history of the 
LMC, the production of the warp in the MW disk and the 
formation of the MS are discussed in §5.41 We conclude 

in m 

2. EARLIER MODELS 

The present galactocentric 3D space, tangential and 
radial velocities, positions, past orbital period, perigalac- 
ticon and apogalacticon distances for the LMC as deter- 
mined by MF80, GSF94, GN96, HR94, vdM02 and M05 
are listed in Table [TJ These works serve as characteristic 
examples of the range of orbital parameters previously 
considered for the Magellanic system. 

With the exception of M05, the listed authors all 
modeled the MW as a singular isothermal sphere and 
accounted for the effects o f dynamical friction using 
the Chandrasekhar formula (jBinnev fc Tremaind (1987) 
equation 7-18) with a constant Coulomb logarithm of InA 
in the range 1 — 3. M05 modele d the dark matter h alo of 
the MW using an NFW profile (jNavarro et al.lll997f l with 
concentration parameter c = 1 1 , virial radius R V i r = 200 
kpc and a virial mass of M vir = 1.1 x 10 12 M Q . However, 
they also considered a substantially lower 3D velocity for 
the LMC (250 km/s). The listed authors were able to 
find solutions where the Clouds follow quasi-periodic or- 
bits that slowly decay owing to dynamical friction/tidal 
stripping as the Clouds move through the dark matter 
halo of the MW. 

The 3D space velocities and lcr errors for the LMC 



as determined by Kl are listed in row four of Table [TJ 
K2 follow the MF80 prescription to search for solutions 
within their proper motion error distribution in which 
the LMC and SMC maintain a stable binary system for 
roughly a Hubble time: a characteristic set of these solu- 
tions is listed in row five of Table[TJ(see Figure 12 in K2). 
Although their analysis illustrates that the new proper 
motion measurements allow for orbital histories that are 
roughly consistent with previous models, the apogalac- 
ticon distances have increased substantially: the LMC 
now reaches distances >150 kpc on its previous orbital 
passage, compared to distances of ~100 kpc determined 
by previous authors. 

Using the mean Kl 3D velocities (row four of TableQJ, 
the apogalacticon distance increases to 220 kpc and the 
orbital period to three Gyr. This contrasts starkly with 
the orbit corresponding to the vdM02 averaged value, 
which yields results similar to those of GN96. The full 
orbital history of the LMC estimated using the Kl mean 
values and those of vdM02 and GN96 are plotted for 
comparison in Figure [TJ (blue solid line, black dotted line 
and red dashed line, respectively). Only the orbital pa- 
rameters of HR94 are consistent with the new picture. 

Although the Kl mean values still imply that the 
LMC has completed several pericentric passages about 
the MW within a Hubble time, the increased orbital pe- 
riod and apogalacticon distance present substantial chal- 
lenges to current tidal and ram pressure stripping models 
for the formation of the Magellanic stream (see £|5.4.3[) . 
Moreover, it is unclear that an isothermal sphere model 
is appropriate at distances >200 kpc. This is certainly 
true over long timescales since dark matter halos evolve 
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Fig. 1. — Galactocentric radial distance plotted as a function 
of time for the mean Kl (solid, blue line), GN96 (dashed, red 
line) and vdM02 (dotted, black line) values. The MW dark mat- 
ter halo is modeled as an isothermal sphere and the effect of dy- 
namical friction is approximated using the Chandrasckhar formula 
with Coulomb logarithm InA = 3. Even in this idealized case, the 
apogalacticon distance reached using the mean Kl values is ~ two 
times larger than the GN96 or vdM02 results and the orbital period 
has increased to ~ three Gyr. 



over time ( Wcchslcr ct al. 20021) For example, the work of 
iWechsler et aiT (l2002) suggests that the formation time 
of a MW-type halo (defined as the time when the halo 
was half as massive as today) is ^8 Gyr ago. If the MW 
halo were smaller and less massive in the past, the or- 
bital energy of a satellite would increase as we int egrate 
its orbit backwards in time (|Peharrubia et ai] [2004) . The 
apogalacticon distance is thus expected to increase more 
dramatically in the past than shown in Figure Q] In light 
of this, we opt for a more detailed MW model as de- 
scribed in the next section. 

3. METHODOLOGY & MODEL PARAMETERS 

Using the current 3D velocity and position of the LMC 
the differential equations of motion {I} can be solved nu- 
merically, allowing the position and velocity of the LMC 
to be followed backwards in time (see MF80). We con- 
sider only the gravitational influence of the MW and 
dynamical friction owing to the passage of the LMC 
through the dark matter halo of the MW: 



d 

r = ^-^mw (|r|) 
or 



' DF 



Mlmc 



(1) 



where Mlmc is the mass of the LMC, r is its position 
vector, 0mw is the potential of the MW and Fdf is the 
dynamical friction term. We trace the orbital history of 
the LMC from equation flU) usi ng the symplectic leap frog 
integration scheme outlined in Spri ngel et al.l (|2001l ). 

Instead of an isothermal sphere model, we describe the 
MW as a static, axi-symmetric, 4-component model, con- 
sisting of a dark matt er (DM) halo, Miyamoto- Nagai disk 
(|Mivamoto fe NagaHll975D , Hernquist bulge (|Hernquistl 



1990) and a halo of hot gas in hydrostatic equilibrium 
with the DM. 



MW 



DM Halo 



HotHalc 



disk 



bulge 



(2) 



The DM halo is initially modeled as an NFW halo. We 
then use CONT RA, a publica l ly ava ilable code written 
by Oleg Gnedin iGnedin et all (|2004f ), to model the adi- 
abatic contraction of the NFW halo in response to the 
slow grow th of an exponent i al dis k. CONTRA follows 
either the Blumcn thal et all (|1986f ) algorithm, which as- 
sumes spherical symmetry, homologous contraction and 
circular particle orbits, or a modified algorithm devel- 
oped by (|Gnedin et alj [20041 ) that attempts to account 
for the high orbital eccentricities of particles predicted in 
hie rarchical structure forma tion scenarios. We use only 
the lBlumenthal et alj (|1986f ) algorithm for simplicity: as 
we demonstrate in ^4.3.11 the model parameters do not 
substantially change our results. The Miyamoto-Nagai 
profile is the monopole expansion of the exponential disk 
profile and is found to be an adequate approximation. 

Our fiducial model is consistent with both model Ai 
of KZS02 (M vir = 1O 12 M ) and known observational 
constraints. Characteristic parameters for the fiducial 
model are summarized in column two of Table [3] The 
third column, labeled High Mass, will be discussed in 
§ 14.3.11 In all cases we adopt a flat ACDM cosmol- 
ogy with parameters consistent with KZS02: h = 0.7, 
&m = 0.3 and baryon fraction f — Q^/C Im = 0.1. Ac- 
cordi ng to the spherical top hat model (|Gunn fc GottJ 
119721) . the virial radius (R V i r ) is defined as the radius 
where the average mass enclosed (pm) equals the virial 
density (p V i r = A^pm)- We follow the analysis of 
KZS02 and take A„j r 340 for our cosmological model. 
Assuming a static MW halo, R V i r is defined as: 



Rvir = 206/i 1 kpc 



M 



97.2 



-1/3 / M . x 1/3 



10 12 h^M, 



( 3 ) 

The resulting rotation curve within 20 kpc is plotted 
in Figure [2 The blue horizontal lines indicate the 3D 
velocity estimates of M05, GN96 and Kl. Notice that 
the M05 velocity is close to the circular velocity whereas 
that of Kl is ~ the escape velocity at 50 kpc; i.e. at the 
current location of the LMC. Triangles repr esent observa- 
tional constraints from HI measurements bv lKnapp et al.l 
(1985), where the solar radius is assumed to be eight kpc. 
Known observational con straints and predictions from 
the MW mo dels of KZS02,|Dehnen fc Binnevl (jl998l ) and 
Smit h~et alj (|2006l ) are compared to our models in Ta- 
ble [4j all values are shown to be consistent with our 
fiducial model. 

Note that if the MW were treated as a point mass, 
the escape velocity at any radius should be \J~2 times 
the circular velocit y (v C i rc ). However, as emphasized by 
iSmith et alj (|2006[ ) , the fact that estimates of the escape 
velocity at the solar radius are larger than \f2v C irc (see 
Table 0|) supports the existence of a massive dark matter 
halo. Thus it is also consistent that v esc » 380 km/s 
> \f2v C i rc — 280 km/s at 50 kpc: v esc probes the total 
mass rather than the enclosed mass. 

We describe the gravitational drag induced by the 
perturbed density field resulting from the LMC's mo- 
tion through the DM halo of the MW using the Chan- 
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TABLE 3 

Model Parameters 



Parameter 


Fiducial a 


High Mass b 


Virial Mass M vir (M ) 


10 12 


2 x 10 12 


Virial Radius R v i r (kpc) 


258 


323 


Halo concentration c 


12 


9 


Disk Mass M di8k (M ) 


5.5 x 10 1(1 


5.5 x 10 10 


Disk scale length r disk (kpc) 


3.5 


3.5 


Disk scale height z^isk (kpc) 


rdisk/5-0 


rdisk/5-0 


Bulge Mass M bu i ge (M ) 


10 10 


10 10 


Solar distance Rq (kpc) 


8.0 


8.0 


Mass fraction of baryons in the disk 


0.055 


0.0275 



a The Fiducial model is consistent with model Ai of KZS02. b The 
High Mass model is consistent with model A4 of KZS02 and is dis- 
cussed in § 14.3. II 



TABLE 4 
Model Constraints 



Parameter 


Observations 


Fiducial 


High Mass 


Other models 


Baryonic Mass inside R Q (10 10 M Q ) 


> 3.9 [1] 


4.1 


4.1 


4.5 [2] 


Escape Velocity at Rq (km/s) 


550-720 [3] 


549 


698 


544 [4] 


Total Mass inside 50 kpc (lO n M ) 


3.2-5.5(4.0-6.4) [5]; 


3.8 


5.0 


3.7 [2] 










1.8-5.6 [6] 


Hot halo # density at 50 kpc (10 -5 cm -3 ) 


< 50 [7] 


25 


80 


1 [8] - 10 [9] 


Total Mass inside 100 kpc (1O U M ) 




6.0 


9.1 


6.0-6.6 [10]; 5.8 [2] 


Total Mass within R vir (1O 12 M ) 


1.1-2.2(1.5-3.0) [11]; 


1.0 


2.0 


1.0 [2]; 0.88-2.56 [4] 




1.5 [12]; 0.2-5.5 [6] 









References. — [l] |Binnev & ; jEvana l|200lD : constraints from the number of microlensing eve nts within Rq ; [2] KZS02: 
model Ai; [3] Browriet aU d2 006D : estimates from hypervelocity stars; [4 Smith ct al. (2006): quoted is their median 
likelihood; [5] Kochanck (1996): constraints from t he local escape velo c ity of stars, disk rotation curve and the motions 
of satellite galaxies, exclu dingf including) Leo I; [6 1 [Wilkinson Evana l|1999T) : motions of satellite galaxies and globular 
clusters at r > 20 kpc; [7] [Rasmusscn & Pcdcrscn upper limit for the ce ntral halo density of the Local Group; [8] 

|Murali|j2000T) : lower limit from the lifetime of the MS; [9] [M oorc & Davis (199j|) : estimated from ram p ressure arguments 
for the formation of the MS; [10] Dchncn &c Binncy ( 1998): their models 1-4; [11] Sakamoto ct al. (2003): constraints from 
kinem atics of MW satellites, globular clusters and horizontal-branch stars, excluding(including) Leo I; [12] Dchncn ct al. 
l |2006|) : determined from radial velocity dispersion profiles of halo objects. 



drasekhar formula J?]): 



DF = 



4^G 2 MLln(A)p(r) 



erf(X) 



2X 



=exp(-X 2 ) 



v 
(4) 

where p(r) is the density of the host halo at the Galac- 
tocentric distance of a satellite of mass M sat , v is the 
orbital velocity of the satellite and X = v/\^2a. Here, 
a is the ID velocity dispersion of the DM halo, which 
can be determined from the Jeans equation; we adopt 
the analytic approxim a tion fo r an NFW profile derived 
by IZentner fc Bullockl (|2003l ) (their equation 6). Al- 
though the Chandrasekhar formula is strictly appropri- 
ate for the idealized scenario of an infinite, uniform, 
non self-gravitating stellar background, dynamical fric- 
tion studies suggest that this app roximation can be 
appli e d to more general sit uations (Ze ntner fc Bullockl 
20031 iTavlor fc Babull Most previous studies of 

the orbital evolution of the Clouds have described the 
Coulomb logarithm as a c onstant value A = &max/&min 
(|Binnev fc Tremainel Il987f ). where 6 max is a fixed cut- 
off radius and & m i n is the im pact parameter. I nstead , 
we adopt the prescription of Has himoto et al.1 (2003), 
wherein & max is replaced by the radial position of the or- 
biting satellite at any given time and 6 m i n = lAk, where 
k = 3 kpc is the softening length if the LMC were mod- 
eled using a Plummer profile. We show in the next sec- 



tion that, contrary to the results of all previous models, 
the influence of dynamical friction on the orbital evolu- 
tion of the LMC is minimal for our fiducial model since 
the LMC has only recently entered the high density re- 
gions of the MW halo. Consequently, uncertainties in 
the analytic prescription for Fdf do not affect our con- 
clusions. Also, throughout this analysis we ignore the 
effects of tidal stripping since this process is similarly 
inefficient at large distances from the halo center. 

4. PLAUSIBLE ORBITAL HISTORIES 

In the following sections we present the results of our 
orbital analysis for our fiducial MW model and discuss 
possible constraints from the position of the MS on the 
plane of the sky. We further illustrate that our results 
are robust to changes in model parameters: namely, the 
mass and shape of the MW halo. 

4.1. Discussion of Fiducial Model Orbits 

Here we present plausible orbital histories for the LMC 
following the methodology outlined in the previous sec- 
tion. In our analysis, 10,000 combinations of the west 
and north components of the proper motion (/j-w^n) 
were randomly sampled within the error space of Kl's 
measurements. For each of these combinations, the 3D 
velocity was derived and the orbit of the LMC was fol- 
lowed backwards in time for a Hubble time using our 
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Fig. 2. — The rotational velocity profile for our fiducial MW 
model is shown within 20 kpc. The circular velocity profile for each 
model component is also plotted individually. Horizontal blue lines 
indicate the 3D velocity estimates of GN96, M05 and Kl. The Kl 
value is close to the red line indicating the escape velocity at 50 
kpc. T he triangles represent the HI measurements of[Knapp ct al. 
(1985) and were used to constrain the rotation curve within the 
solar circle (8 kpc). 



fiducial MW model. All other parameter errors, e.g., 
that of the distance modulus and the current position, 
are small in comparison to those of the proper motions 
and do not affect this analysis. 

In Figure [3] the orbital path of the LMC is traced back- 
wards over four Gyr and plotted in the Galactocentric YZ 
plane. The evolution of the galactocentric radius of the 
corresponding orbits are plotted as a function of time in 
Figured! where t = corresponds to today. In both fig- 
ures, the bold-face dashed(dotted) line indicates the orbit 
traced by the old GN96(vdM02) velocities in our fiducial 
model. Comparing these lines to their counterparts in 
Figure [TJ it is immediately apparent that the orbital his- 
tories of the Clouds in a cosmologically-motivated MW 
model are dramatically different from those in an isother- 
mal sphere model. In our fiducial model, the vdM02 
and GN96 proper motions (Tabled]) imply that the LMC 
has completed only one orbit within 10 Gyr and reached 
an apogalacticon distance of 300-400 kpc, whereas these 
same proper motions suggest orbital periods of 1.5 Gyr 
and an apogalacticon distance of 100 kpc in the isother- 
mal sphere model. These striking results indicate that, 
independent of the proper motion measurements, the 
choice of MW model can significantly alter our picture of 
the orbital history of the Magellanic system. The depen- 
dence of these results on model parameters are discussed 
in $Lp] 

In Figure[3]we also plot the orbital path of the LMC for 
each of the 10,000 proper motion combinations allowed 
within Kl's error space. We find that the orbits are still 
roughly polar, but are not strictly confined within the 
YZ plane. Colored areas indicate the portion of the YZ 
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Fig. 3. — Orbital paths corresponding to 10,000 randomly se- 
lected combinations ((j-w^n) traced four Gyr in the past and plot- 
ted in the Galactocentric YZ plane. The present location of the 
LMC is indicated by the black square and the MW disk plane by 
the short-dashed line. The dash-dotted line indicates the virial 
radius of the halo (258 kpc). Orbital paths corresponding to fiw 
values within ±lcr from the mean value (p^ = —2.03) are confined 
within the red area. Paths that correspond to fiw components 
within (1-2)<t are confined within the light blue region, between 
(2-3)tT within the blue region and between (3-4)<r within the green 
region. Only orbits within the upper green region (nw = t^w +4(t) 
cross the disk plane within 4 Gyr. The bold-face dashed and dot- 
ted lines indicate the orbits traced by the GN96 and vdM02 values, 
respectively, in our fiducial model. 



plane spanned by orbits corresponding to (hw,Hn) com- 
binations within a specified a from the mean (/i*). The 
orbit corresponding to the mean values (/i* = (-2.03,0.44)) 
is parabolic. Orbits with a larger \fiw\ component (-la 
to -4cr from /j,^) correspond to hyperbolic orbits. The 
magnitude of | must decrease by 3-4er before the orbit 
will cross the disk plane. In Figure [4] these same orbits 
are constrained between the red and green lines, which 
correspond to orbits with fi = fi* ± 4c. Even in the case 
most favorable to bound orbits (/j,w — t^w + ^c), the 
orbital period is roughly a Hubble time and reaches an 
apogalacticon distance of 550 kpc, which is substantially 
larger than that predicted by all other studies. On such 
distance scales dynamical friction plays little role in mod- 
ifying the orbital history of the LMC and uncertainties 
in the analytic description of Fdf are irrelevant (§ [3]). 

We now examine these results more statistically: in 
Figures [5] and [6] all 10,000 combinations are color coded 
according to the time (Figure [5]) or distance (Figure O 
at which the LMC last crossed the MW disk plane. In 
both cases the light blue dots indicate orbits that never 
crossed the MW disk plane within tn- Dashed lines indi- 
cate the number of standard deviations of a given point 
from the mean proper motion ((/xjy,^^); black triangle). 
The large asterisk shows the vdM02 average of previous 
proper measurements: this value is well outside 4er of 
li^y, but is consistent with [i* N . The distance and time 
of disk crossing for the vdM02 value (250 kpc; 1.5 Gyr 
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Fig. 4. — Orbital evolution of the LMC plotted as a function 
of time in the past, where Time=0 corresponds to today. The 
dash-dotted line indicates the virial radius of the halo (258 kpc). 
All allowed orbits are constrained by the red and green lines, which 
indicate the outer boundary of the green regions in Figure^^^ = 
li^v ± 4<t). The blue line represents the orbital evolution of the 
LMC if the mean values are used (/i* = (—2.03,0.44)). Even in 
the best case scenario, the LMC completes only one orbit within 
tn and reaches an apogalacticon distance of 550 kpc. The dashed 
and dotted lines indicate the orbits traced by the old GN96 and 
vdM02 values, respectively, in our fiducial model. 

ago) are marked in the legend. 

No solutions within la of the Kl mean ever cross the 
disk plane: in most of those cases the LMC never com- 
pletes a single orbit within tn ■ The 2% of cases that do 
cross the disk plane do so at substantially larger times 
and distances than predicted by any previous study. In 
fact, only ~0.1% of all cases ever cross the disk plane 
< 4 Gyr ago; i.e. at times greater than twice the orbital 
period predicted by GN96 using an isothermal sphere 
MW model. These plots thus provide a measure for the 
timescale and strength of the interaction between the 
MW and LMC. 

From Figures [5] and [6j the dominant factor controlling 
the orbital history of the LMC is the west component of 
the proper motion (fj,w) : if & solution exists for a given 
fxw value, then it will exist for all jiijv. The west compo- 
nent determines the tangential velocity - as the magni- 
tude of nw increases, so does v tan and, correspondingly, 
the orbital eccentricity. As such, the results for proper 
motion estimates with /ijv values beyond the scale of the 
y-axis can be extrapolated from the result of any proper 
motion combination with the same iiyj ■ For example, the 
GN96 proper motion estimate of (hw,Hn) — (-1.72,0.12) 
is well outside 4er of Kl's mean values and does not ap- 
pear in our figures. But, by the above argument, we 
expect the GN96 results to fall somewhere between the 
+4<t results and those of the vdM02 proper motions. 
Reading from Figures [5] and if \iw « —1.72, the GN96 
orbit should have crossed the disk plane ~ 1.7—3 Gyr 
ago at a distance of ~ 260 — 400 kpc. Checking Figures [3] 
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Fig. 5. — 10,000 points randomly sampled from the (4<r) proper 
motion error space of the Kl measurements for the LMC. The 
dashed ellipses indicate the standard deviation of the enclosed 
points from the mean (black triangle). For each point, the orbital 
history of the LMC was computed by integrating the equations of 
motion backward in time for the fiducial MW model. Solutions 
that cross the MW disk plane at least once within tjj are color 
coded by the time of the last disk plane crossing. Light blue dots 
indicate solutions that never crossed the disk plane within 
There are no solutions that cross <3 Gyr ago and in only ~0.1% 
of all cases does the LMC cross the disk plane < 4 Gyr ago; i.e. 
twice the orbital period predicted by previous studies. The lack of 
colored points less than fiyy ~ —1.9 indicates that the west com- 
ponent of the proper motion is the dominant factor in determining 
the orbital path of the LMC. The asterisk indicates the average of 
previous proper motion measurements as determined by vdM02 - 
the corresponding orbit crossed the disk plane 1.7 Gyr ago. 

and 21 we find that indeed the GN96 orbit crossed the 
disk plane at a distance of ~300 kpc, ~2 Gyr ago (the 
GN96 orbit spans only the YZ plane). 

In Figure [7] the colored (disk crossing) solutions from 
the previous figures are re-plotted as a function of the 
time and radius of crossing. They are now color-coded 
by the standard deviation of the proper motion from the 
mean value. Even within 4ct of the mean, the LMC only 
crosses the disk plane at timescales longer than three Gyr 
and at distances larger than 415 kpc. Moreover, it is 
unreasonable to consider the LMC as an isolated system 
or ignore cosmology over such distances and timescales. 

We thus conclude that in our fiducial model the or- 
bit of the LMC must be close to parabolic: the LMC is 
on its first passage about the MW and is currently at 
perigalacticon. 

4.2. Orbital Constraints from the MS 

We now ask whether we can rule out either our fiducial 
(NFW) or isothermal sphere MW models by comparing 
the projection of the corresponding orbits of the LMC on 
the sky to the current posi tion of the MS. This ex ercise 
is motivated by the work of lJohnston et all (|1999bf l. who 
suggest that tidal streams act as "fossil records" of the re- 
cent orbital history of their progenitors and could thereby 
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Fig. 6. — Same as Figure [5] except that solutions that cross the 
MW disk plane at least once within tjj (only 2% of all cases) are 
color coded by the radial distance at which the LMC last crossed 
the disk plane. Light blue dots again indicate solutions that never 
crossed the disk plane. There are no solutions that cross the disk 
plane at distances smaller than 400 kpc; i.e., at four times the 
apogalacticon distance predicted by previous studies. The orbit 
corresponding to the vdM02 values crossed the disk plane at a 
distance of 260 kpc. 

provide a probe of the galactic potential. We note that 
the formation mechanism of the MS is still a subject of 
debate; however, both the ram pressure and tidal strip- 
ping models predict that the orbit of the Clouds (approx- 
imated here just by the LMC) should trace the MS for at 
least some time in the past (e.g. GN96 and M05). Since 
the isothermal sphere and the fiducial (NFW) model give 
very different orbital histories in 3D space, it is worth ex- 
ploring whether they give drastically different projected 
locations for the LMC in comparison to the MS. 

Figure [8] shows the LMC orbits corresponding to the 
Kl and GN96 velocities for both our fiducial and isother- 
mal models as seen by an observer situated at the loca- 
tion of the Sun looking towards the South Galactic Pole. 
The orbits are traced backwards in time until the LMC's 
line of sight position vector subtends an angle of f00° 
with respect to its present-day posi tion vector: this cor- 
responds to the extent of the MS (|Putman et alj 120031 . 
hereafter P03). 

GN96 specifically chose 3D velocities for the LMC such 
that its orbit in an isothermal sphere MW model would 
trace the observed position of the MS (solid green line). 
We find little deviation in the projected orbital path of 
the LMC if our fiducial model is used with the GN96 ve- 
locities (dashed green line) , even though the orbits differ 
substantially in 3D space (see Figures Q] and |4j . 

The dashed (solid) red line traces the LMC's orbit us- 
ing the mean velocities of Kl for the fiducial(isothermal 
sphere) model. Once again, in projection there is little 
deviation between the two MW models. There is, how- 
ever, a striking difference between the allowed Kl orbits 
and the GN96 results: the LMC no longer traces the MS 
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Fig. 7. — This plot summarizes the results of Figures \5\ and [6] 
Here, all solutions that crossed the disk plane (the colored points 
from the previous two figures) are plotted as a function of the 
time and radial distance at crossing. Points are color coded by 
the deviation of the proper motion combination from the mean 
value: /i refers to some combination of (J-w an d £tjv> although in 
practice only fiyy matters. Only 0.1% of cases within l-2a of the 
mean ever cross the disk plane within tn . Also, no orbits cross the 
disk plane < 3 Gyr ago or at distances smaller than 415 kpc. The 
symbol ©(*) indicates the result for the GN96(vdM02) values in 
the fiducial model. 



(solid green line). 

The dashed (solid) blue lines represent the worst and 
best fits to the position of the MS for the fidu- 
cial(isothermal sphere) model. Orbits allowed within 4ct 
of the mean are bounded by those lines, as indicated by 
the black arrow. These bounds were determined by ex- 
tremizing the x 2 residuals of the LMC's orbit to a linear 
parametrization of the P03 HI data: Magellanic latitude 
B = 0.04L + 4.0. The Magellanic coordinate system 
(L,B), defined such that its equator is parallel to the MS 
and the center of the LMC is currently located at Lw- 
40 degrees, has been introduced he re as a simple way 
to vi sualize the Magellanic system (|Wannier fc Wrixonl 
[1971). 

The good agreement between the projected orbits for 
our fiducial and isothermal sphere models imply that we 
cannot use the location of the MS to distinguish between 
MW models. Although a comparison of the linc-of-sight 
velocities of the orbit versus the HI data for the MS could 
serve as a better discriminant, the correlation is strongly 
dependent on the formation mechanism of the MS (see 

EH). 

We have ignored the SMC so far in this analysis, 
even though previous models have assumed that the MS 
actually represents material torn from the SMC (e.g. 
GN96). The purpose of this section is two- fold: firstly, to 
show that we cannot differentiate between the NFW and 
isothermal halo models based on the current position of 
the MS; and secondly, to illustrate that the proper mo- 
tion measurements necessitate a revision of current the- 
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Fig. 8. — Possible orbits of the LMC are mapped as a polar 
projection in galactic (l,b) coordinates. The orbit is followed back- 
wards in time from the LMC's current position (black triangle) 
until it extends 100° in the sky. Solid(dashed) lines indicate orbits 
computed using the isothermal sphere(fiducial) model. The green 
lines trace the orbit implied by the GN96 velocities: GN96 chose 
their initial velocities such that the LMC's orbit traces the the true 
position of the stream (P03). The red lines show the values using 
the mean Kl velocities, while the blue lines depict the best and 
worst fits of the allowed orbits to the position of the MS (max and 
min x 2 ! indicated by the black arrow). 

ories for the origin of the MS. The inclusion of the SMC, 
located 20 degrees from the LMC on the plane of the sky 
(or a separation of 23 kpc), does not affect either of these 
conclusions. 

As we discuss in H5.ll the SMC's impact on the or- 
bital history of the LMC is minimal, owing to the large 
LMC:SMC mass ratio of -10:1. Thus, the global dynam- 
ics of the LMC alone is sufficient to determine whether 
this analysis can differentiate between halo models (see 
also j)4.3p . Furthermore, if we require that the Clouds 
have been bound to each other for at least some time in 
the past, the orbital path of the SMC will be strongly 
dependent on that of the LMC. Consequently, a shift in 
the LMC's orbit with respect to the current location of 
the MS will manifest as a shift in the orbital path of the 
entire Magellanic system. As such, the deviation cannot 
be explained away by assuming that the MS originated 
instead from the SMC (see Figure [H] in HB~T|) . 

We may also worry that this pronounced deviation im- 
plies an error in the Kl measurements. As mentioned 
in <JTJ HR94 consider a similarly large tangential veloc- 
ity component as Kl, yet their fiducial orbit traces the 
MS. This is because, like GSF94 and GN96, they a priori 
chose /zat to ensure alignment. But these ^jv estimates 
are substantially different from not only the new Kl val- 
ues, but also the weighted average of previous measure- 
ments as determined by vdM02 (see Table [5]). Specifi- 
cally, if we ignore the Kl values and repeat this analy- 
sis using the vdM02 proper motions instead, we get the 
same result: the vdM02 estimate of ^jv is similar to that 



of Kl. Figure [9l shows the projected orbital path of the 
LMC using the Kl and vdM02 proper motions (red and 
blue lines, respectively) overlaying P03's HI data for the 
MS. The Kl and vdM02 results are nearly indistinguish- 
able, illustrating that the clear deviation in the projected 
orbital path and current location of the MS cannot be 
dismissed as an error in the new measurements. 

As such, regardless of the choice of MW model, there 
are severe implications for either the tidal or ram pressure 
stripping theories for the origin of the MS: the LMC's or- 
bit using the Kl mean velocities deviates from the cur- 
rent location of the MS by ~7 degrees on the sky, or ^7 
kpc if the MS traces a circle of radius 55 kpc. This devi- 
ation might be attributable to a number of second order 
effects independent of the formation mechanism of the 
MS, e.g, the rotation of the MW halo or the L/SMC's 
disk might affect the direction in which the material is 
removed. Regardless, such effects do not change our bot- 
tom line: all previous theoretical studies of the orbital 
evolution of the Magellanic Clouds assumed that the MS 
and the LMC's orbit are co- located on the sky, mean- 
ing they are inconsistent with both the new HST and the 
weighted average of all previous proper motion measure- 
ments. The reliability of the Kl measurements will be 
discussed in detail in H5.2I 

4.3. Model Dependences 

Here we show that our analysis is robust to changes in 
model parameters. Specifically, we examine whether the 
conclusions of £ 14.11 change if we were to consider a MW 
model with either a non-spherical or more massive halo 
than in our fiducial case. 

4.3.1. Mass Dependence 

We again consider a static 4-component MW model 
and maintain the same cosmological parameters. We also 
keep the same disk and bulge parameters, but increase 
the virial mass to 2 x 10 12 M Q . The halo concentration 
is decreased to c = 9 in ord e r to m atch the HI obser- 
vational data of iKnapp et al.l (|1985f ). Model parameters 
are summarized in Table [3] in the column labeled "High 
Mass" and the resulting rotation curve is plotted in Fig- 
ure |T0j This model is consistent with KZS02's model 
A 4 (M vir = 2 x 1O 12 M , c = 10). It also satisfies the 
observational constraints as listed in Table IH although 
the escape velocity at th e solar radius is higher than the 
estimates of lSmith et ail (|2006l ) (498 < v esc < 608 km/s) 
and the local number den sity of the gaseous halo c o nflicts 
with the upper limits of iRasmussen fc Pedersenl (|200lD 
(see Table |4]). Note also that the escape velocity at 50 
kpc is now 552 km/s, which is substantially higher than 
the mean velocity of Kl. 

Higher mass models require even lower concentrations 
so that the maximal circular velocity does not exceed 
the observed value. For example, if the upper limit of 
M vir — 5.5 x 10 12 M (|Kochaneklll996[ ) is adopted, a con- 
centration of c = 6 — 7 is required. KZS02 rejected higher 
mass models based on the assumption that the concen- 
tration parameter should be in the range c = 10 — 17; if 
M V i r > 2 x 1O 12 M and c > 10, the maximum circular 
velocity would violate the HI observations. This range for 
the concentration parameter follows from ACDM model 
predictions that the virial mass a nd the concentratio n 
parameter are strongly correlated ijBullock et alj [2001). 
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Fig. 9.— The HI distribution of the MS from the data of P03 is 
plotted as a polar projection using the same scale as that in Fig- 
ure [8] The axis are the galactic longitude (gl) and latitude (gb). 
The SMC is located 20 degrees away from the LMC, corresponding 
to a physical separation of 23 kpc. The green line traces the orbit 
theoretically determined by GN96 in an isothermal MW model. 
Note that since GN96 advocate for a scenario in which the MS 
originates from the SMC , the green line does not match the data 
perfectly (see Figure Il9l for a plot that includes the SMC's orbit). 
The blue(red) line shows the path traced by the vdM02(Kl mean) 
proper motions for the fiducial MW model. The vdM02 result is in- 
distinguishable from that of the Kl mean values and both deviate 
by ~ 7° from the location of the MS: the observationally deter- 
mined proper motions differ markedly from theoretical estimates. 

Moreover, for the iKochanekl (|1996l ) upper mass limit, 
the escape velocity at the solar radius is ~ 1000 km/s 
- in comparison, the fastest hypervelocity star observed 
to date has a Gala ctic rest-frame velocity of 709 km/s 
(jBrown et al J 12006( 1. 

We repeat the analysis of Sj4] using our High Mass 
model: the orbits of the LMC implied by the same 10,000 
randomly selected fj,w an d combinations are traced 
backwards for a Hubble time. The Galactocentric radial 
position of the LMC is plotted as a function of time in 
the past in Figure 1111 which complements Figure 3] for 
the fiducial model. Although the results have improved, 
the orbital period for the mean values (/i*; blue line) is 
~6 Gyr and the LMC reaches an apogalacticon distance 
of ~400 kpc. Thus, even in this High Mass model, the 
LMC travels on a highly eccentric orbit and is unlikely 
to have undergone more than one pericentric passage. In 
the best case (n^y + 4cr; red line in Figure ITTj) the LMC 
completes three orbits within the past 10 Gyr, but the 
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Fig. 10. — Rotation curve plotted within 20 kpc for our High 
Mass MW model (M vir = 2 X 10 10 M ). The bulge and disk com- 
ponents are the same as in the fiducial model, while the halo mass 
has been increased by a factor of two. The NFW halo concen- 
tration parameter was decreased relative to the fiducial model in 
order to match the HI observational data within the solar radius 
( Kna pp et a CUD; triangles). The escape velocity is now 552 
km/s at 50 kpc, which is higher than the mean velocity estimate 
of Kl (top blue line). This model is consistent with KZS02's max- 
imum halo mass model A4. 



apogalacticon distance is still > 150 kpc. The orbit corre- 
sponding to the GN96 velocities in this High Mass model 
is shown for comparison (dashed black line) . Both these 
cases are approaching the isothermal sphere solutions of 
GN96, shown as the dashed red line in Figure [TJ 

We also repeat the statistical analysis of $2J Figures H"2l 
and [13] complement Figures [5] and El respectively. As 
before, orbits that cross the disk plane are color coded 
by the time or distance of the last crossing. As expected, 
there are more solutions that cross the disk plane at least 
once within % (80% of cases cross the disk plane <4 
Gyr ago). But, solutions still cannot be found for large 
\fiw\ components (black dots) and solutions within la 
of the mean do not cross the disk plane < 2.4 Gyr ago 
or at distances <300 kpc, which corroborates our earlier 
results. 

We conclude that, if the LMC is moving with a veloc- 
ity close to the mean value determined by Kl (v = 378 
km/s) and if the MW is not well approximated as an 
isothermal sphere at distances > 200 kpc, the LMC is 
unlikely to have undergone more than one pericentric 
passages about the MW within a Hubble time. 

4.3.2. Axial Ratio of the MW halo 

So far, we have considered only spherical halo models: 
if we vary the axial ratio, can a better fit between the 
position of the MS on the sky and our projected orbits 
be obtained? Given the good correlation between the 
projected orbits for our fiducial and isothermal sphere 
models, we assess the effects of halo sphericity using the 
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Fig. 11. — Same as Figure [4] except that the orbital evolution is 
traced using the High Mass model. The dash-dotted line indicates 
the virial radius of the halo model (323 kpc). All allowed orbital 
histories are constrained by the red and green lines, which corre- 
spond to /i^y ± Act. The shortest allowed orbital period is ~ 3 Gyr 
(red line). The orbit corresponding to the mean proper motions 
(blue) has an orbital period of ~6 Gyr and reaches an apogalacti- 
con distance of 400 kpc on the last passage. The orbit using the 
GN96 velocities in this model is also plotted for comparison (black 
dashed line); it is approaching the isothermal sphere solution of 
GN96 (see Figure Q). 

isothermal sphere model of Binncy (1981): 



* = ^log(i* 2 + * 2 /V) 



(5) 



where we have adopted cylindrical polar coordinates 
(R,4>,z). The parameter q defines the axial ratio of 
the halo potential; in general, the corresponding density 
profile will be more non-spherical. We vary q between 
q = 0.5 (oblate) and q = 1.5 (prolate; spherical halos cor- 
respond to q = 1) and compute the resulting orbits for 
the mean Kl velocities 1 . The orbits are plotted in pro- 
jection in Figure [TH We find that prolate halos provide 
better x 2 fits than oblate halos. However, the position of 
the orbit relative to the MS remains largely unaffected 
by the halo shape in the recent past - which is where the 
tidal and ram pressure stripping models predict the best 
agreement. 

Th ere are few const r aints on the sphericity of the MW 
halo. iFellhauer et al.l (|2006| ) determined that the angu- 
lar differences on the sky between the bifurcated stellar 
tid al streams of the Sagitt arius dwarf spheroidal detected 
bv lBelokurov et aD (|2006| ) implied minimal precession of 
the orbital plane. They concluded that the axial ratio 
of the halo must be close to spherical. However, this 
model does not explain the apparent velocity dif ference 
of ~ 15 km/s between the bifurcated streams: iHelmil 
( 2004) suggests that the stream kinematics favor prolate 

1 We have checked the corresponding density profiles and ensured 
that they are reasonable. The oblate halos have the most non- 
physical density profiles and are slightly dimpled at the poles. 
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Fig. 12. — Same as Figure [5] except that orbits are computed for 
the High Mass model. Now 80% of orbits cross the disk plane <4 
Gyr ago: most orbits cross the disk plane at least once within tjj, 
although if /ny < (flyy ~ 3<r) solutions still cannot be found (black 
dots). Orbits within ltr do not cross the disk plane < 2.4 Gyr ago 
and no orbits cross the disk <1 Gyr ago. Using the vdM02 values, 
the LMC crosses the disk plane 0.81 Gyr ago. 

halos, as determined from radial velocities of M giants in 
the leading arm selected from 2 MASS (Two Mic r on Al l 
Sky Survey). On the other hand. [Johnston et al.l (2005) 
used the 2MASS M giants to trace the tidal streams and 
determined that the level of precession was best matched 
by oblate halos. Clearly there are no definitive answers. 
Unfortunately, Figure [T3] shows that our analysis cannot 
provide additional constraints on the halo sphericity. 

We go one step further and explore models in which, 
in addition to varying the axial ratio, we also vary the 
slope, j3, and amplitude, 77, of the rotation curve at ~ 50 
k pc. We use as our basis the "power law" galaxy model 
of lEvand (|1994h in which the potential is given by 



(R 2 



(6) 



where k 



p [ i? LMC+ z LMc/9]' i/2 - This prescription 
gives a series of uninterrupted galaxy models with rising 
(ft < 0) and falling (/3 > 0) rotation curves at large 
radii and properties abutting those of the isothermal 
sphere. We explore models with —0.7 < (3 < 0.7 and 
0.9i>o < 77 < 1.1 «o- However, after an extensive param- 
eter search, neither of the parameters (3 or 77 improve 
the agreement with the MS. The parameter q is found to 
make the biggest difference to the shape of the orbit. 

Since model parameters do not appear to be able to 
explain the deviation between the projected position of 
the MS and the orbits determined by the Kl velocities, 
a revision of current theories for the formation of the MS 
may be warranted. The implications of our analysis to 
formation mechanisms of the MS is continued in detail 
in 35X51 

5. DISCUSSION 
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Fig. 13. — Same as Figure [5] except that orbits are computed 
for the High Mass model. Orbits within lcr do not cross the disk 
plane at Galactoccntric distances < 300 kpc and no orbits cross 
at radii < 150 kpc. Using the vdM02 values, the LMC crosses the 
disk plane at a distance of 120 kpc. 



In the following sections we show that our conclusions 
are robust to perturbations from the SMC ( §5.1|) and con- 
sider caveats to our analysis: namely the reliability of the 
recent proper motion measurements ( ^5.2p . We further 
discuss the likelihood of a first passage scenario ( §5.3|) 
and the implications of the new velocity measurements 
for the star formation history of the LMC, the nature of 
the warp in the MW disk and the origin of the Magel- 
lanic stream (E i5.4j) . Specifically, we illustrate that even 
if the MW were modeled as an isothermal sphere, there 
are significant challenges to our current understanding of 
the formation of the Magellanic stream (MS). 

5.1. The SMC-LMC Binary System 

The LMC:SMC mass ratio of -10:1 likely precludes 
the SMC from being a major determinant in the LMC's 
orbital evolution. However, observations of a common 
envelope of HI gas surround the L/SMC and a bridge of 
material linking the two Clouds imply that the L/SMC 
must have maintained a binary state for some time in 
the past. In this section we discuss the role of LMC- 
SMC interactions to the robustness of the predicted LMC 
orbits and the long-term stability of a binary LMC-SMC 
system. 

We again follow the methodology outlined in ££3j how- 
ever the equations of motion for the LMC have been 
modified from equation |T]) to include the gravitational 
interactions between the SMC and LMC. We model the 
LMC (SMC) using a Plummer profile with a softening pa- 
rameter of 3(2) kpc. The SMC mass is assumed to be 
2 x 1O 9 M (SMCTMC mass ratio of 1:10). The equa- 
tion of motion for the SMC is analogous, except that we 
also include a dynamical friction term that acts on the 
SMC when it enters the current tidal r adius of the LMC 
(-15 kpc; Ivan der Marel et all (|2002f) ). We follow the 
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Fig. 14. — Possible orbits of the LMC using only the mean Kl 
values mapped as a polar projection in galactic (l,b) coordinates 
for isothermal sphere MW models with varying axial ratios. The 
orbit is followed backwards in time from the LMC's current posi- 
tion (black triangle) until it extends 100° in the sky. The green 
line traces the orbit implied by the GN96 velocities in an isother- 
mal model of the MW. The other lines represent orbits in halos of 
varying sphericity: light blue (q=0.5; oblate), red (q=1.0; spheri- 
cal) and blue (q=1.5; prolate). 



analysis of IBekki fc Chibal (l2005f) and app roximate this 
friction term as ( Binnev fc Tremaindll987l ): 



F LS = -0.428 ln(A 



GM§ v LS 



LS 



VLS 



(7) 



where i>ls i s the relative velocity between the Clouds and 
the Coulomb logarithm is ln(AL,s) = 0.2. 

Here, we repeat the statistical analysis of S]4] with the 
inclusion of the SMC. 10,000 pairs of points were drawn 
randomly from both the LMC and SMC proper motion 
error space (K1,K2). The corresponding SMC and LMC 
orbits were traced backwards in time for a Hubble time 
in our fiducial MW model. In Figures [TS] and [T|5] we 
re-plot Figures [5] and [7] now including the SMC. The 
SMC introduces more scatter into the plots and there 
are now more solutions within 2er, but they exist at large 
distances and over long timescales. Thus, the overall 
results are not substantially changed: there are still no 
orbits within lcr that cross the disk plane. 

Although there are large uncertainties in the 
LMC:SMC mass ratio, increasing the SMC mass by a 
factor of a few does not substantially change these re- 
sults. The same is true if this analysis were repeated 
using the High Mass mo del described in £14.3.1 

The calculations of IBekki fe Chibal (|2005l) for an 
isothermal sphere MW model suggest that it is unlikely 
that the LMC and SMC remained a bound system for a 
Hubble time since dynamical friction owing to the SMC's 
passage within the LMC's halo d estab ilizes the binary 
system. In support, IBekki et al.l (|2004f) argue that the 
different cluster formation histories of the LMC and SMC 
may be explained if they originated as separate entities 
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Fig. 15. — Same as Figure[5]except that interactions between the 
LMC and SMC have been included. The SMC introduces scatter, 
but the previous results are unchanged: there are no solutions that 
cross the MW disk plane within la of the mean or for large l/ivH- 



(see also £15 .4|) . However, given the sparse distribution 
of the outer MW satellites, a capture scenario of the 
SMC by the LMC seems improbable. We now exam- 
ine whether the LMC and SMC can maintain a stable 
binary state in our fiducial MW model if the LMC is 
moving at the Kl mean velocity of 378 km/s. 

10,000 points were randomly drawn from the SMC 
proper motion error circle (K2) and plotted in Figure [T7l 
they are color coded based on the longevity of the binary 
state. Bound orbits are identified by whether the SMC's 
velocity relative to the LMC (wls) is smaller than the 
local escape velocity. In most cases this criterion is not 
initially satisfied (light blue dots; t<0.1 Gyr); i.e. the 
SMC is currently not bound to the LMC. But solutions 
where the binary state is maintained for 8 < t < <h 
(red squares) are allowed within lcr of the mean. Thus, 
although dynamical friction between the clouds was in- 
cluded, binary states that last for a Hubbl e time are plau- 
sible, which is contrary to the results of iBekki fc Chiba 
(2005). This is because the relative velocity between the 
LMC and SMC is determined by the tidal force imparted 
by the MW, and our fiducial model diff ers significantly 
from t he isothermal sphere model used bv lBekki fc Chiba 
(2005). Correspondingly, very few stable binary states 
exist when this analysis is repeated for our High Mass 
model and it becomes easier to maintain the binary state 
if the LMC: SMC mass ratio is increased. This is consis- 
tent with a similar analysis by K2 in an isothermal sphere 
model (i.e. a high mass model), in which they found < 
10% of orbits to be in a stable binary state. 

Repeating this analysis using the limiting 4cr LMC ve- 
locities (corresponding to fj,w — ji^ ± 4<r) instead of the 
mean LMC velocity, we find that a stable binary state 
(tbound > 8 Gyr) can always be located. However, it 
becomes increasingly difficult to identify stable binary 
states as the \fiw\ component is increased. The require- 
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Fig. 16. — Same as Figure [7] but including interactions between 
the LMC and SMC. The plot is more scattered but still consistent 
with the previous results: there are no cases that crossed the disk 
plane <2.3 Gyr ago or at distances <350 kpc and the best results 
are obtained if fi > (fi* +3<r). The symbol ©(*) indicates the result 
for the GN96(vdM02) values. 



ment of binarity could thus be used to further constrain 
the SMC or LMC proper motion error space: in order 
for the L/SMC to form a binary system their transverse 
motions must be comparable, and the new transverse 
motions for the LMC have increased substantially. In 
Figure [T8l the orbital paths of a stable binary system for 
the mean and limiting LMC velocities are traced over a 
Hubble time in the Galactocentric YZ plane. We con- 
clude that for every LMC orbit that is allowed by the 
data, there also exists a bound SMC orbit that is also 
allowed by the data. 

Considering only long-lived binary orbits (red squares 
in Figure [IT)) , we now repeat the analysis of £|4.2l and con- 
firm the statement that the inclusion of the SMC cannot 
explain the discussed deviation in the orbital path of the 
Magellanic system with respect to the current location of 
the MS. In Figure [T9l the projected path followed by the 
SMC (dashed line) and LMC (solid line) using the old 
GN96 values in an isothermal sphere model (green) are 
plotted in Galactic coordinates. Notice that the SMC's 
path traces the true location of the MS (shown in Fig- 
ure [9|) better than the LMC does. However, for the al- 
lowed Kl LMC velocities, the GN96 estimate of the SMC 
proper motion (indicated by the symbol 0) does not yield 
a stable binary LMC-SMC system. Instead, we plot the 
most stable binary state if the LMC is moving at the 
mean Kl values in the isothermal sphere (red) and fidu- 
cial (blue) MW models. Here the fiducial model orbits 
also correspond to the orbits traced by the blue lines in 
Figure [T71 As expected, the deviation of the LMC's or- 
bit translates into a shift in the orbital path of the entire 
Magellanic system. Thus, if the LMC and SMC formed 
a binary system for some time in the past, the LMC's 
orbit provides a fair assessment of the behavior of the 
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Fig. 17. — 10,000 points randomly drawn from the SMC proper 
motion error ellipse (K2). Each point is color coded by the amount 
of time the SMC remains bound to the LMC, i.e. how long 
Vesc > fLSi if the LMC is moving at the Kl mean velocity. In 
the majority of cases, this criterion is not currently satisfied (light 
blue dots). The red squares indicate solutions where the binary- 
system is maintained for tpj > t > 8 Gyr: such solutions do exist 
within la of the mean (black triangle). The symbol Q indicates the 
proper motion corresponding to the GN96 estimate of the SMC's 
current velocity: in our model, this value corresponds to a currently 
unbound LMC-SMC system. 



entire Magellanic system. 

5.2. Reliability of the New Proper Motions 

There are three new ingredients in this analysis: Kl's 
fJ-w, l^N and the MW mass model. The far fewer number 
of perigalactic passages for the Clouds is a generic pre- 
diction of our fiducial model for the MW and is essen- 
tially independent of the proper motion data. However 
the conclusion that the Clouds arc on their first pas- 
sage relies on the Kl fi\y value and corresponding error 
bar. Given the unexpected implications for the history 
of the MW-LMC-SMC system (see jgg) , it is worth re- 
examining the new (j,yy value in the context of other work. 
Kl's value for fiw of the LMC is on the high end of the 
values reported previously (see Figure 14 in Kl), but it 
is consistent within the e rrors with the Hippar cos i value 
(jKroupa k. Bastianlfl997h and the (jPedreros et al.ll2002ft 
value; see also Table [2] 

The line-of-sight velocities of LMC carbon stars can 
also be used to put constraints on the transverse veloc- 
ity of the LMC center of mass. Even though the rota- 
tion curve and the transverse velocity component are not 
uniquely determined by the line-of-sight velocity field, 
plausibility arguments based on the Tully-Fisher relation 
and comparison to M33 puts the transverse component 
between 200 and 600 km/s (vdM02). Thus the transverse 
velocities used in this analysis are not ruled out by the 
carbon star kinematics. An increase in the east motion 
of the LMC gives better agreement with the line-of-sight 
velocity field than the vdM02 proper motions, as can be 
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Fig. 18. — Stable binary orbits plotted in the galactocentric 
YZ plane for the mean (blue) and limiting LMC velocities, which 
correspond to fiw = A 1 ^ ± 4<r (red and green). The orbits are 
traced over a Hubble time. The LMC orbit is indicated by the 
solid line and the SMC by the dashed line. The orbital path of the 
LMC is not significantly affected by the SMC and a stable binary 
state can be found within the SMC proper motion error space for 
all allowed LMC velocities. 



seen from Figure 8 in vdM02. The dotted trapezoid in 
this figure is the 68.3% confidence region obtained from 
the line-of-sight velocity field under the assumption that 
the rate of change of inclination of the LMC disk, di/dt, 
is zero. In Kl a combined analysis of the carbon star 
kinematics and the proper motion field is used to show 
formally that di/dt is more consistent with zero than 
with the use of previous (lower) proper motion values. 

In a recent analysis comparing the kinematics of 
HI, carbon stars and red supergiants in the LMC 
(01 sen fc Massevl (|2007l) ). the authors were able to cor- 
rect the line-of-sight velocities for the Kl proper motion 
and make a useful interpretation of the resulting veloc- 
ity fields, finding that each of these tracers defines a flat 
rotation curve. The large tangential component of the 
LMC velocity does not introduce spurious motions in the 
analysis of their line-of-sight data. Finally, a substan- 
tial tangential component for the LMC was predicted by 
iLin &: Lvnden-Belll (|1982l ). in which they explore a tidal 
model for the origin of the MS in a massive dark halo 
model for the MW. Tangential velocities consistent with 
these values were also predicted by HR94 in a ram pres- 
sure stripping model for the MS (also with a massive 
MW halo). 

Notably, Kl's estimate of hn is consistent with that 
of vdM02 (see Table [2|) . This implies that the apparent 
shift in the orbital path of the Magellanic system with 
respect to the current location of the MS, as discussed in 
£14.21 exists even if the old proper motion measurements 
are used (see Figure [9]). Theoretical models have invoked 
smaller values of /i^ in order to satisfy the underlying as- 
sumption that the Clouds should trace the MS. However, 
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Fig. 19.— Possible orbits of the LMC (solid lines) and SMC 
(dashed lines) are mapped as a polar projection in galactic (l,b) 
coordinates. The orbits are followed backwards in time from the 
Clouds' current positions (black triangles) until they extend 100° 
in the sky. The green lines trace the orbits implied by the GN96 ve- 
locities in an isothermal model of the MW (Iso). The red and blue 
lines show the projected orbits for the most long-lived binary state 
if the LMC is moving at the mean Kl velocities in the isothermal 
and fiducial MW models, respectively. 

this assumption is strongly dependent on the formation 
mechanism of the MS, which in itself is highly uncertain 
(see ggX3]) . 

5.3. Evidence for a First Passage Scenario 

The possibility that the LMC is on its first pas- 
sage ab out the MW is not a novel i dea. As 
noted bvlFich fc Tremamel (|1991hlLin et alj (| 19951) and 
ISawa & Fuiimotol (|2005) , if the LMC is currently at peri- 
galacticon, its orbital angular momentum is comparable 
to the rotational angular momentum of the MW's stellar 
disk. Since the LMC's orbit is near polar, the spin axis of 
the MW disk is perpendicular to that of the LMC's orbit. 
Thus, if the LMC has been in a decaying quasi-periodic 
orbit about the MW for a Hubble time, the origin of 
its high angular momentum is puzzling. The new higher 
tangential velocity measured by Kl only exacerbates this 
"angular momentum problem" and cannot be explained 
b y current models of the orbital ev olutio n of the Clouds. 

iRavchaudhurv fc Lvnden-Belll (|1989f ) suggest that 
tidal torques exerted by M31 early in the LMC's or- 
bital history may be responsible for the L MC's hi g h tan - 
gen tial velocity . Buil ding on this theory, IShuterl (|1992l ) 
and lBvrd et alj <| 1994f) considered the possibility that the 
Clouds underwent a close encounter with M31 and were 
only recently tidally captured by the MW. Although our 
computed orbits do not support a scenario where the 
LMC originates from M31, the source of the LMC's high 
tangential velocity may be easier to understand if the 
LMC was not always bound to the MW. In this way, 
the evolution of the LMC-SMC-MW system cannot be 
considered in isolation, and the LMC-SMC binary may 



have been subjected to external torques before entering 
the MW's virial radius. 

In addition, based on a morphological comparison of 
the sa tellite galaxies of the MW and M31. lvan den Berghl 
(|2006l ) suggests that the L/SMC may be interlopers that 
originally formed in a more remote region of the Local 
Group. All inner satellites of M31 are early-type objects, 
which is also true of the MW except for the presence of 
the irregular LMC and SMC at small galactocentric radii 
with res pect to the other MW s atellites. Moreover, as 
noted bv lFich fc Tremaind (|1991f ). the max imum line of 
sight v elocity of the MS (wlsr = -410 km/s: lBruns et all 
(2005)) is not only higher than that of the Clouds them- 
selves (262 ± 3.4 km/s; vdM02) but also of any satellite 
within 200 kpc. Thus, if the orbital path and velocities 
of the LMC are similar to that of the MS, they conclude 
that its orbits must be very different from those of other 
satellites (see also the disc ussion in 8 15 .4.3|) . 

Following the analysis of lLin et al.l ( 19951 ) we next con- 



sider simple consistency checks for a scenario where the 
LMC is currently unbound to the MW. If the LMC's or- 
bital eccentricity is negligible, the dynamical mass of the 
MW within the orbit of the LMC can be approximated 
as follows: 
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(8) 



where vtan is the tangential velocity of the LMC (367 
km/s; Kl) and -Rlmc = 50 kpc is the current Galacto- 
centric distance of the LMC. This is much larger than 
observational constraints (see Table [4]) and so the LMC 
cannot be on a circular orbit. 

Since the radial velocity is non-zero (v ra( i = 89 km/s; 
Kl), we know the LMC must be on an eccentric orbit. 
Assuming that all the galactic mass is within 50 kpc, a 
lower bound on the dynamical mass can be estimated as: 
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where r a is the apog a lactic on distance. Following the 
analysis of iLin et al.l (l995), if the LMC is bound to 
the MW then the apogalacticon distance must be less 
than the Galactic tidal radius with respect to M31, 
so r a <300 kpc. The corresponding mass estimate is 
roughly the upper limit on the mass estimates within 50 
kpc (6 .6 x 10 11 M Q within 50 kpc; Wi lkinson &: Evani 
(1999)). This simple lower bound suggests that the LMC 
is likely on a highly eccentric orbit and could be effec- 
tively unbound. 

A first passage scenario for an LMC-type galaxy at 
the present epoch is not at odds with current theories 
of the hierarchical build-up of dark matter (DM) halos. 
In a study using both N-body simulations and Extended 
Press-Schechter calculations to estimate the fraction of 
MW size halos that have experienced recent mergers, it 
is found that 70% have accreted an LMC-sized object in 
the past —5 Gyr (James Bullock, private communication 
2007). Here "accreted" means "first falls within" a -300 
kpc virial radius. For our fiducial model, the LMC typi- 
cally enters the virial radius within the past 1-2 Gyr. The 
corresponding fraction of such mergers in the simulations 
is 30%. We thus conclude that a first passage scenario 
for an LMC-type object is not statistically improbable. 
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Finally, it should be noted that the primary constraint 
on the long-term orbital history of the Clouds comes from 
current formation mechanisms of the MS, which require 
that the Clouds have undergone multiple pericentric pas- 
sages, and from theoretical links between the observed 
star formation rates and close encounters between the 
LMC-SMC-MW system. These issues are explored in 
detail in the next section. 

5.4. Phenomenological Implications 

In the following we consider the implications of the 
new proper motion measurements and consequent orbital 
history of the Clouds for our understanding of their star 
formation histories, the formation of the warp in the HI 
layer of the MW disk, and the origin of the MS. 

5.4.1. Star Formation History of the Magellanic Clouds 

The star formation history (SFH) of the Magellanic 
Clouds is believed to be strongly influenced by dynami- 
cal an d hydrodynamical galaxy interactions (Wcstcrlund 
1997). Detailed knowledge of the orbital history of the 
Clouds is thus critical to the understanding of their 
SFHs. In the picture we have developed in §|4] (our fidu- 
cial MW model), the Magellanic clouds are on their first 
passage about the MW and entered within the virial ra- 
dius of the MW DM halo only - (1 - 3) Gyr ago. As 
they travel supersonically (mach ~3) through the halo 
gas, they interact hydrodynamically with the ambient 
medium, forming a bow sh ock that will increase the gas 
density at the leading edge (|de Boer et al.|[l998f h As the 
Clouds approach the galactic center they are also sub- 
ject to tidal forces from the MW. We therefore expect 
the Clouds to exhibit a heightened SFR within the past 
few Gyr. 

The SFH of the LMC's disk is believed to be relatively 
smooth and continuous, whereas t he recent SFH of its 
bar is more episodic. Specifically, ISmecker-Hane et all 
(2002) identified the dominant stellar populations in the 
bar with episodes of star formation that occurred 4-6 
and 1-2 Gyr ago. These episodes have been attributed 
to tidal interactions during pericentric passages between 
the LMC, SMC and the MW. Moreover, an "age gap" 
has been observed in the LMC's globular clusters, which 
refers to the identification of only o ne cluster in the 
LMC with an a ge between ~3-13 Gvr (TRich et alj|2001t 
iDa Costal Il99l . iBekki fc Chibal (|2005j) note that the 
reactivation of cluster formation ~3 Gyr ago could be 
explained by the onset of hydrodynamic interactions be- 
tween the LMC's disk and the ambient halo gas. The 
resulting bow-shoc k could act as a cat alyst for cluster 
formation (see also lde Boer et all (|1998f l). If the Clouds 
have undergone multiple pericentric passages about the 
MW, it is unclear why such hydrodynamical effects would 
be important only recently - whereas this picture is a 
natural consequence of a first passage scenario. 

A rise in the SMC's m ean SFR within the past 3 Gyr 
has also been inferred by [Harris fc Zaritskvi 1)20041 ) . who 
describe the recent star formation rate as highly irreg- 
ular with bursts of star formation occurring at 2.5, 0.4, 
and 0.06 Gyr ago. Unlike the L MC, the SMC is un - 
barred and has no spiral structure (|Zaritskv et al.ll2000f ). 
thus the mechanism for driving star formation in the 
SMC likely involves galaxy interactions. Specifically, 



IZaritskv fc Harris] (|2004D suggest that the irregular ap- 
pearance of the SMC principally owes to recent star for- 
mation triggered by hydrodynamic interactions between 
the SMC's gaseous components and the halo gas (or outer 
gas envelope of the LMC), rather than tidal effects. Thus 
the recent SFH of the SMC also appears to support a first 
passage scenario. 

5.4.2. Warp of the MW Disk 

From 21-cm observations it is clear that the HI layer 
in the MW disk is substantially warped (vertically dis- 
torted to distan c es >3 kpc above the galactic plane; 
Diplas & Savage (1991)). The dust layer has also been 
observed to be similarly displaced from the galactic plane 
(|Freudenreich et ai1ll994l ). Although the origin of the 
warp is still unknown, the LMC has historically been 
considered as providing a probable excitation mecha- 
nism owing to its proximity and mass. However, even 
if the LMC were on a circular orbit of radius 50 kpc, 
its tidal fields are not strong e nough to distort the HI 
MW disk by > 2 kpc ve rtically (|Burkdll957t lKerrlll957t 
iHunter fc Toomrel I1969D . This could be improved if 
the LMC were substantially more massive than previ- 
ously estimated (> 10 n M Q ; iBinnevi (|1992h ) or if the 
pericenter of its o rbit w ere much smaller (~20 kpc; 
IHunter fc Toomrel (|1969t )1. Our analysis suggests that 
the LMC is on a parabolic orbit and is only now at its 
closest approach to the MW: on such an orbit, the torque 
exerted by the LMC is insufficient to excite a warp of the 
observed magnitude. Thus, if the Magellanic Clouds are 
the source of the warp, their influence on the MW needs 
t o be amplif i ed. 

IWeinberel (|1998f ) proposed that the warp results from 
a combination of two processes: 1) direct tidal forcing by 
the LMC; 2) resonant forcing of disk bending modes by 
the density disturbance excited by the LMC's passage 
through the halo. If the LMC fo llows a non-deca ying 
quasi-perio dic orbit, as assumed bvlWeinberj (|1998l ) and 
recently by IWeinberg fc Blitz] (|2006l ). it will excite dis- 
crete modes in the halo. The higher order modes de- 
cay quickly, whereas low order modes persist over a few 
dynamical times; moreover, on a quasi-periodic orbit 
the weakly damped modes will be continu a lly re- excited 
and so dominate the response. iTsuchival (|2002t ) tested 
this scenario by means of N-body simulations and de- 
termined that if the LMC follows a decaying orbit, it 
could excite a warp of the required magnitude over a 
6 Gyr timescale. Note t hat the N-body simulations of 
iGarcia-Ruiz et all (|2002t ) yielded only a 25% amplifica- 
tion of the total torque exerted by the LMC, rather than 
the r equired factor of five to reproduce the observed 
warp. IWeinberg fc Blit j (|2006h suggest that the discrep- 
ancy owed to an unluck y set of parameter choices by 
IGarcia-Ruiz et al.l (|2002T ). 

Altering the satellite's orbit changes the forcing fre- 
quency, and thereby affects both th e amplitude and 
the orientation of the warp response. ITsuchival ((20021 ) 
adopted a MW model and orbital parameters such that 
the LMC's orbital period is ~1.5 Gyr and the maxi- 
mal apogalacticon distance is ~ 100 — 110 kpc. But, 
for an isothermal sphere model, the Kl mean velocities 
imply an orbital period of ~ 3 Gyr. The LMC thus 
completes only 2 orbits wit hin 6 Gyr (Figu re [1]) , rather 
than 4-4.5 orbits as in the ITsuchival (|2002f ) model, and 
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reaches a maximum apog alacticon distance of 250 kpc. If 
iGarcia-Ruiz et all (|2002l ) found null warps owing to their 
halo parameter choices, the drastically different orbital 
values suggested by the new observations will surely limit 
the effectiveness of this mechanism. 

In the picture we have developed here, the LMC may 
be on its first passage about the MW. Fly-by encoun- 
ter s have been considered in th e context of halo wakes 
by iVesperini fc Weinberg) {2000). Here, the interloping 
galaxy induces a continuous spectrum of frequencies in 
the halo. Since only the low order modes are long-lived, 
the net effect is similar to that produced by a satellite 
on a quasi-periodic orbit. However, the magnitude of 
the effect is strongly dependent on the impact parameter 
and speed of the encounter: the disturbance is significant 
only after the perturbing galaxy has reached pericenter 
and the ma gnitude decreases as its veloc ity is increased. 
Specifically. IVesperini fc Weinbe rg (2000) suggest that a 
perturber that is 5% of the primary's mass and mov- 
ing at ~200 km/s could substantially perturb the pri- 
mary's halo during a fly-by encounter. Since the new 
measurements imply that the LMC (mass ~2% of the 
MW) is traveling at ~ 400km/s and we have concluded 
that the LMC is only currently at pericenter, it seems 
unlikely that the density perturbations induced by the 
LMC's passage through the MW halo can provide suffi- 
cient torque to excite a warp of the observed magnitude 
in the MW disk. 

5.4.3. Formation of the Magellanic Stream 

The origin of the MS is highly controversial and dif- 
ferent studies advocate for a variety of formation scenar- 
ios. Historically, the MS has been thought of as a tidal 
feature, where gas is stripped from either the LMC or 
SMC at their last pericentric passage. But since tidal 
forces should work equally on both stars and gas, the 
lack of any stellar tidal feature associated with the MS 
(jGuhathakurta fc~R citzcl 1998) suggests that it may not 
have a purely tidal origin. Many works have also con- 
sidered the role of hydrodynamic processes, such as ram 
pressure stripping. However, both of these mechanisms 
rely on the assumption that the L/SMC have undergone 
multiple pericentric passages about the MW. 

In the tidal stripping scenario, a sate llite loses mass 
predominantly at per i centric passages dJohnston et all 
119951: Uohnstonl 119981 : Uohnston et~atl Il999aft In the 
GN96 picture, the MS formed from material stripped 
from the SMC at its previous perigalactic approach 1.5 
Gyr ago, which also coincided with a close encounter with 
the LMC. Our computed orbits (Sj4j do not support this 
scenario: for the isothermal sphere model the previous 
pericentric passage occurred ~ 3 Gyr ago; for our fidu- 
cial models we have determined that the L/SMC are cur- 
rently at pericenter and only in the "best case scenario" 
(/■*tV+4<t) wm another pericentric passage occur within a 
Hubble time. Without multiple pericentric passages, the 
effectiveness of tidal stripping is severely limited; along 
the fi ducial orbits, the instantaneous tidal radius ( King 
fl962l) of the LMC(SMC) is greater than 15(7) kpc un- 
til roughly 0.2 Gyr ago, when the Clouds were ~70 kpc 
away from the Galactic center. But the bulk of the HI 
in the LMC is distributed within a radius of ~5 kpc 
(|Stavelev-Smith et alJl2003D and the SMC has a line of 
sight depth of ~7 kpc (|Harris fc Zaritsk~vll2004h . In ad- 



dition, material outside the tidal radius is not strippe d 
instantaneously (|Johnstonl fl998: Joh nston et al.| [T999a). 
It is therefore doubtful that a tidal stripping scenario can 
produce the MS if the Clouds have only been subjected 
to significant tidal effects over the past 0.2 Gyr. Fur- 
thermore, the deviation of the projected orbits and the 
current position of the MS in the recent past presents a 
challenge to tidal stripping models, as discussed in ^4.21 
During high velocity encounters, hydrodynamic inter- 
actions with the ambient medium are likely more sig- 
nificant than tidal distortions. Specifically, various ram 
pressure stripping scenarios have been proposed to ex- 
plain the formation of the MS, e.g.. iMoore fc David 
(|1994 hereafter MD94), HR94 and M05. Owing to 
the low density of the ambi ent MW halo gas (n e = 
10~ 5 - 10~ 4 , see Table H the iGunn fc Gottl (fl972T ) for- 
malism for ram pressure stripping is insufficient to re- 
move material instantaneously from the LMC disk. This 
is true even if the LMC's disk were face-on (which is not 
the case - the LMC disk is inclined ^30 degrees from 
the line of motion) or if Kl's new larger velocities were 
adopted. 

Alternatively, MD94 proposed that the MS was ram- 
pressure-stripped from the Clouds when they last passed 
through the extended ionized disk of the MW. They sug- 
gest this collision occurred ~500 Myr ago at a galactocen- 
tric radius of 65 kpc. At this distance the density in the 
disk may be high enough for ram pressure to overcome 
the gravitational restoring force in the SMC or intracloud 
region. However, this scenario is not supported by our 
orbital analysis: for our computed orbits the LMC either 
never crosses the disk or does so at large radii (>400 kpc) 
and at much earlier times (>3 Gyr ago) (see Figures [5] 
and[6|). Even for our high mass model the orbits do not 
support a disk crossing at radii <150 kpc (see Figure ITS"]). 

Along a different line, M05 proposed a "continuous 
stripping" scenario, wherein the LMC gas disk edge is 
heated as a result of the compression exerted by the 
ambient halo gas. They argue that the corresponding 
increase in thermal energy can unbind the gas in the 
LMC disk on time scales comparable to the orbital time. 
This model implicitly assumes that the LMC has under- 
gone multiple pericentric passages, which M05 achieved 
by using a significantly lower velocity than had been pre- 
viously estimated (250 km/s). For our orbits, the inter- 
action time between the LMC and the halo gas (< 1 — 2 
Gyr) is likely too short for this to be a viable mechanism. 
Other hydrodynamic effects, such as Kelvin-Helmholtz 
and Rayleigh- Taylor instabilities (e.g. viscous stripping, 
iNulsenl {1982)), operate over even longer timescales and 
are expected to contribute minimally. 

Alternate scenarios have been proposed where loosely 
bound material stirred by on-going tidal interactions 
within the LMC-SMC binary system is stripped from 
the Clouds or intracloud region (e.g. HR94 and 
lYoshizawa fc Noguchil (|2003f )). Such mechanisms may be 
viable, as they are less dependent on the assumption of 
multiple pericentric passages. We have also identified 
long-lived binary LMC-SMC states that are allowed by 
the data for all our computed orbits (see However, 
such scenarios will still have difficulties explaining why 
the Clouds' orbits do not trace the MS (E2J). 

Our analysis suggests that the origin of the Magel- 
lanic stream remains an open question: all previously 
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proposed mechanisms need to be revisited in light of the 
surprising consequences of the new proper motion mea- 
surements of Kl a nd K2. 

As suggested by iFich fc Tremainel (|1991[ ) , the orbital 
history of the Clouds should be similar to that of the 
MS. Ultimately, a way to distinguish between forma- 
tion mechanisms may thus come from the ability of the 
proposed mechanism to account for the deviation of the 
LMC's orbit with respect to the location of the MS on 
the sky (Figure [5]). Furthermore, it is also expected that 
the line of sight orbital velocities with respect to LSR 
(ulsr) should be roughly consistent with that of the HI 
observed in the stream today. 

In Figure [201 the wlsr is computed at each point 
along the orbits depicted in Figure [8] and plotted as 
a function of Magellanic longitude (L). The black line 
is the HI velocity data of P03, estimated from the ze- 
roth moment (intensity, Int) and first moment (veloc- 
ity) radio maps. At each Magellanic longitude (L), the 
average over the Magellanic latitude (B) is plotted as 
< Mos >= J2 v ios{B)Int(B)dB/J2lnt{B)dB. Once 
again, the green line indicates the results for the isother- 
mal GN96 model, which approximates the HI velocity 
data. 

The i>lsr values computed along the orbit correspond- 
ing to the mean Kl velocities in both the isothermal 
sphere (solid red line) and fiducial models (dashed red 
line; allowed values are bounded by the blue lines as in- 
dicated by the black arrow) are substantially higher than 
the observations of P03. The tidal stripping model of 
IConnors et alj (120051) is the only w ork to date to account 
for the recent iBruns et all (|2005t ) HI PARKES survey 
results; however, their orbital velocities are inconsistent 
with the new proper motions. The details of the pre- 
dicted velocity structure of the MS are strongly depen- 
dent on the formation mechanism and further discussion 
is beyond the scope of this analysis. This will be the 
subject of a forthcoming paper. 

6. CONCLUSIONS 

This study has been motivated by two considerations: 
1) previous works regarding the orbital evolution of 
the Magellanic Clouds have considered mainly isother- 
mal sphere models for the MW. However, over the past 
decade our understanding of the structure and formation 
of dark matter halos has become more detailed and we 
can consider perhaps more realistic models to describe 
objects located at large Galactocentric radii. 2) The- 
oretical models for the MS have historically been used 
to constrain the proper motions of the Clouds - with 
the improved observations, we can now assume that the 
proper motions are known to some degree of accuracy. 
These two factors have very specific and surprising im- 
plications for our understanding of the orbital evolution 
of the Clouds. 

We have shown that for a ACDM-motivated model of 
the MW that is consistent with observational constraints 
and most theoretical models (such as model Ai of KZS02; 
i.e., our fiducial model, M vir = 10 12 M©), proper motion 
estimates prior to the recent HST measurements result 
in drastically different orbital histories than previously 
predicted. Specifically, the GN96 estimate of the LMC's 
proper motion (based on theoretical models of the lo- 
cation and line of sight velocities of the MS) and the 
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Fig. 20. — Line of sight velocities with respect to LSR (i>lsr) 
plotted as a function of Magellanic longitude (L) along the orbit 
shown in Figure [8] The lines are color coded as in Figure [8] The 
black line indicates the HI data of the MS from P03 (extending 
until MSV as denned by P03 in their figure 5. The i>lsr magni- 
tudes from our orbit calculations are much higher than the P03 HI 
observations or the GN96 results (dashed green line). The black 
triangle indicates the current systemic LMC line of sight velocity 
(262.2 ± 3.4 km/s; vdM02). 

weighted average of pre-2002 ground-based and Hippar- 
cos proper motion measurements (vdM02), both yield 
orbits where the LMC has completed only one passage 
about the MW within 10 Gyr and reached apogalacticon 
distances of 300-400 kpc, compared to the orbital period 
of 1.5 Gyr and apogalacticon distance of 100 kpc advo- 
cated by GN96 in an isothermal model. Thus, a drastic 
revision of our understanding of the orbital history of the 
Clouds is warranted in such a ACDM motivated picture 
for the MW. This result is independent of the proper 
motion measurements. 

The new data give us additional information: the com- 
ponents of the proper motion control very specific orbital 
parameters of the LMC. The west component determines 
the tangential velocity and thereby sets the number of 
pericentric passages (orbital period), apogalacticon, lo- 
cation of disk crossing and the stability of the binary 
system. The north component dictates the location of 
the orbits with respect to the current position of the MS 
on the sky. The implications of the new data for the Mag- 
ellanic system are discussed below, with specific reference 
to the role of the individual proper motion components. 

The HST proper motion measurements (Kl) yield 3D 
velocities that are substantially higher (~100 km/s) than 
previously estimated, owing to the increase in the west 
component. This implies that the LMC is currently 
traveling at the local escape velocity. We subsequently 
searched the error space within 4<r of the mean values and 
found that the LMC could complete at most one orbit 
within a Hubble time, and only in cases where the west 
component of the proper motion was small (|/x*y| — 4cr). 
We tested whether our results were robust to changes 
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in model parameters by considering a higher halo mass 
model (M v ir — 2 x 10 12 M©) that was still consistent 
with observational constraints. Although the number of 
pericentric passages increased, the orbital period implied 
by the mean velocities is still longer than five Gyr. Per- 
turbative effects owing to the presence of the SMC were 
shown to have a negligible effect on our orbital analysis 
owing to the LMC:SMC mass ratio of ~10:1. Further- 
more, our computed orbits do not preclude the existence 
of a stable binary LMC-SMC system, however this be- 
comes more difficult to maintain as the MW mass in- 
creases. 

We conclude that either the Magellanic Clouds are 
on their first passage about the MW or that the MW 
DM halo is well-modeled as an isothermal sphere to dis- 
tances > 200 kpc (i.e. substantially more massive than 
2 x 10 12 M Q , although note that this may conflict with 
known observational constraints). Even if the rotation 
curve of the MW is flat out to distances of 200 kpc, the 
apogalacticon distance and period of the LMC's orbit im- 
plied by Kl's mean velocity will have increased substan- 
tially compared to previous theoretical estimates (220 
kpc and 3 Gyr versus 100 kpc and 1.5-2 Gyr). This is a 
direct result of the substantial increase in the /iw proper 
motion component measured by Kl. 

The north component of the proper motion does not 
affect the above result, but as mentioned, it controls the 
location of the orbit when projected on the plane of the 
sky. Using the mean Kl values, the LMC's orbit devi- 
ates from the position of the MS by 7° (~ 7 kpc) on the 
sky. This result is independent of both the (i\y compo- 
nent and the halo model (see Figure [8]). This is signifi- 
cant because while there are theoretical models that have 
predicted tangential velocities as high as Kl, they have 
all assumed /in ~ (e.g., HR94 and LL82) - but this is 
not reconcilable with the data. Even if we ignore the Kl 
measurement (/xjv = 0.44 ± 0.05), the vdM02 weighted 
average of previous measurements (/xjv = 0.34 ± 0.16) is 
not consistent with zero. This presents an impediment 
to both tidal and ram pressure stripping models of the 
MS, since both assume that the MS is co- located with 
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the past orbit of the Clouds. 

In conclusion, if a ACDM-motivated model accurately 
describes the MW, and if the Kl results are correct, 
the Clouds are on their first passage about the MW. A 
first passage scenario is consistent with morphological 
studies of the satellite populations of the MW and M31, 
where the dlrr L/SMC appear as interlopers compared 
to the dSph galaxies that dominate the satellite popu- 
lation at small Galactocentric distances (jvan den Berghl 
2006), and is not statistically improbable from simula- 
tions of hierarchical structure formation. The origin of 
the LMC's high orbital angular momentum may also be 
better understood if the Clouds were not always orbiting 
about the MW. 

A first passage scenario has a number of unexpected 
phcnomenological implications that warrant further in- 
vestigation. It supports a significant increase in the 
star formation rate of the Clouds during the past 1-3 
Gyr, corresponding to when the Clouds first enter within 
the virial radius and begin to interact with the halo 
gas. It also drastically limits the timescale over which 
the Clouds and the MW can have interacted, making 
it highly unlikely that the LMC could have excited the 
warp in the MW disk. This forces a major reassess- 
ment of proposed formation mechanisms for the Magel- 
lanic stream, which depend strongly on the assumption 
that the Clouds have undergone multiple pericentric pas- 
sages. The implications for the formation of the MS are 
severe even if an isothermal sphere model of the MW is 
adopted: the orbital period and apogalacticon distances 
have increased, limiting the effectiveness of tidal and ram 
pressure stripping. 
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